in H «oi Ainio y< $ aoiad Odo 


H65-21667 

- * ' ' " \ 

! 

1 

(ACCESSION 

^7 _ , 

iTHriW# 

7 i 


3Q ; 

—NASA CR OR TMX OR AO NUUuLm 

(CATEGORY) 


X-650-R5-86 

//AsA 



PERIODIC WIND 
COMPONENTS AT 
METEOR HEIGHTS 





BY 

R. G. ROPER 
W. G. ELFORD 



NASA 


MARCH 1965 


GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 



PERIODIC WIND COMPONENTS 
AT METEOR HEIGHTS 


by 

R. G. Roper 
W. G. Elford 


March 1965 


GODDARD SPACE FLIGHT CENTER 
Greenbelt, Maryland 


PERIODIC WIND COMPONENTS AT METEOR HEIGHTS 


R. G. 


by 

Roper* and W. G. Elford 


t 


ABSTRACT 


3 \ui 

The measurement of winds in the height range from 80 to 100 km by means 
of radio reflections from the ionized trails of meteors is a well established 
technique. Until recently, however, the full potential of this technique had not 
been realized. The purpose of this report is to demonstrate the power of an 
analysis developed by Groves (1959), which has been applied to the datacollected 
at Adelaide (35° S) for the thirteen months from December, 1960 to December, 
1961, inclusive. The results show that the most significant periodic zonal and 
meridional wind components are the 24 hour and 12 hour, with a smaller 8 hour 
component occasionally present. In general, the vertical component is random, 
with an amplitude from one to two orders of magnitude less than that of the 
horizontal wind. 
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PERIODIC WIND COMPONENTS AT METEOR HEIGHTS 

by 

R. G. Roper and W. G. Elford 


INTRODUCTION 

The purpose of this report is two-fold, presenting as it does a summary 
of upper atmosphere wind measurements carried out at Adelaide (35° S) during 
1961, as well as documenting the computer programs used in the computation 
and analysis of these winds. The measurements were obtained from the radio 
observation of drifting meteor trails, and have revealed a detailed picture of 
the diurnal and seasonal motion of the atmosphere throughout the height range 
from 80 to 100 km. 

The Meteor Data 

The Adelaide Meteor Project employs a combination of continuous wave 
and pulse techniques to measure the position and line of sight drift of each me- 
teor trail observed. A description of the method has been given by Robertson, 
Liddy and Elford (1953). 

The parameters measured by the system are 

a. the line of sight drift of the trail. 

b. the direction cosines of the reflection point on the trail 
relative to NS and EW axes. 

c. the radar range of the echo point, and 

d. the time of occurrence of the echo. 

The echo rate is not constant with time; an average of 5 useable wind ech- 
oes per hour was obtained during the 1961 survey. To obtain enough data for 
analysis, records were taken for an average of nine consecutive days about 
the middle of each month. The mean height of all echoes for each month was 
93 (±1) km, and reliable winds were obtained for the height range 80 to 100km. 
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The Analysis 

In his meteor wind analysis, Groves (1959) sets up a wind model against 
which the meteor data is matched. The parameters specifying the model are 
the coefficients of the polynomial wind variations with height and time. These 
are conveniently expressed in tabular form. 
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where S is the normalized height given by 

S = (2h = z max - z min )/ (z max - z m iu) 

with z max , z min the maximum and minimum heights respectively over which 
the profile is fitted and h is the height variable within the range z m in, z max* 
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If the variations with height in the prevailing wind and the component of 
period T = 2 it/co are taken as 

a polynomial of order 2 in the EW direction 
" 3 in the NS direction 

and a constant in the vertical (constant with height, not time) 

then 24 parameters (ay, by, cy) are required to specify the motion. 

These parameters are determined by matching the data to the model, using the 
method of least squares. 

Application of the Groves Analysis to the Meteor Data 

No attempt will be made here to present the mathematics of the analysis; 
the reader is referred to Groves (1959). The IBM 7094 FORTRAN IV program 
which performs this analysis of the meteor data is listed as program ERG in 
Appendix I. The number of parameters best specifying the height profiles 
have been arrived at by trial and error. In this context, the parameters 
"best" specifying the profiles are the coefficients of the minimum order poly- 
nomials which adequately represent the height structure of the mean motion. 
Calculating the coefficients of the terms of higher order than necessary is 
extremely wasteful in terms of computer execution time. For example, the 
calculation for a total of 70 periodic components of the coefficients of a cubic 
variation with height of the zonal (EW) component, a quartic variation of the 
meridional (north south), and a constant vertical wind (30 coefficients in all 
for each periodicity) for 1000 line of sight drifts measured over 9 days takes 
42 minutes of execution time on an IBM 7094 Model 2. Assumption of a quad- 
ratic zonal, cubic meriodinal and constant vertical variation (24 coefficients) 
results in an execution time of 22 minutes, a little over half that of the pre- 
vious computation. Attempting to over-specify the variations in the wind can 
also cause the matching process to become unstable, and this has been ob- 
served with (3,4,0) profiles. No such instability has been encountered in any 
of the months so far processed with a (2, 3, 0) profile. 

The Presentation of the Results 


The digital output from program ERG (consisting of the amplitude and 
phase of each periodic component, together with the RMS error in each, for 
each of 16 heights in the range from 75 to 105km) is both listed and written on 
magnetic tape. This tape provides input for the plotter program listed in 
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Appendix II. The program PLOT writes a low density binary tape suitable for 
use on the EAI Data Plotter. All the graphs presented in this report are copies 
of those produced by the EAI machine. 

The following are definitions of the quantities presented graphically in the 
next section. 

a. The wind amplitude spectrum (for a given height). The amplitude of 
each periodic component is plotted against period for each of the three 
wind components (zonal, meridional and vertical). 

b. The "specific" energy spectrum (for a given height). The "specific" 
wind energy is the wind energy per unit mass, calculated for any one 
period r as the sum of the squares of the zonal (u), meridional (v), 
and vertical (w) amplitudes of period r . 

i.e. , 


E t = u T 2 + v T 2 + w T 2 (metres/sec) 2 

Reduction of Actual Data 

An average of 800 echoes for each of the months from December, 1960 to 
December, 1961, inclusive, have been subjected to the analysis outlined in the 
previous sections. 

The following parameters were used as a basis for the analysis for each 
month. 

1. Only those echoes falling within the height range 75 to 105km were 
processed. 

2. Both prevailing and periodic zonal components were allowed up to a 
quadratic variation with height over this height range. 

3. Similarly, both prevailing and periodic meridional components were 
allowed up to a cubic variation. 

4. The vertical component was considered to be constant with height, but 
able to vary with time. 

5. Amplitudes and phases of periodic components were calculated for the 
range of from 0.5 cycles/day to 4 cycles/day (48 hour to 6 hour 
periods) in increments of 0.05 cycles per day. 


4 


The following features have been chosen as being most descriptive of the 
characteristics of the periodic motion in the height range from 80 to 100 km. 

a. The variation of the wind energy spectrum with height, as measured at 
83, 91 and 97km. 

b. The amplitude spectra for zonal, meridional and vertical components 
as measured at the mean monthly echo height of 93 km. 

The six graphs for each month from December, 1960 to December, 1961 
appear as Figures la to 14f, inclusive. 


DISCUSSION 

Insufficient time has as yet been spent on examination of the results to 
justify any detailed interpretation. However, the spectra presented in Figures 
1 to 14 confirm without exception that the dominant periodic components in the 
winds at these levels are the diurnal (24 hour) and semidiurnal (12 hour) com- 
ponents. At the latitude of Adelaide (35° S), the 24 hour is, in general, of 
comparable magnitude to the 12 hour, except at the equinoxes, where the 24 
hour predominates. There is evidence in some months of a significant 8 hour 
component which, when it appears, has an amplitude somewhat less than those 
of the 24 or 12 hour components. 

Although the main concern of this report is to investigate the presence of 
periodic motions, some reference to the prevailing wind components is perti- 
nent. The zonal winds are shown for three heights in Figure 15. The bars 
about each plotted point represent the RMS deviation. It can be seen that over 
the height range 80-100km the zonal wind is predominantly toward the East. 
The only strong wind reversal occurs at the upper level, where the wind is 
toward the West during the winter of 1961. No significant reversals occur at 
the other levels but the wind has its maximum eastwards amplitude in the 
spring at 91km and in the summer at 83 km. As a result of the rapid change 
in the zonal wind with height the seasonal patterns at 83 and 99km are almost 
opposite in phase. This behavior is also reflected in the wind gradients which 
have maximum values of +4 m/sec/km in summer and -4 m/sec/km in winter. 

In contrast to the zonal winds, the meridional winds shown in Figure 16 
exhibit an annual behavior which is similar at all levels. In general, north- 
ward winds occur during summer and southward winds during winter. A sim- 
ilar meridional variation is found for these levels at Jodrell Bank (53° N) and 
Mawson (68°S). Thus the meridional flow at these levels is consistent with a 
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horizontal movement of air from the summer to the winter pole. The Adelaide 
results show that the meridional wind increases with height over the range 80- 
100km and that above 90km the amplitude of the meridional wind is comparable 
to that of the zonal wind. 

It should be remembered, when further interpreting these results, that 
one year's data is not by itself sufficient to justify a detailed analysis on a 
seasonal basis. The results presented here, while considerably more detailed, 
are still complementary to those already published covering the years 1952 to 
1958 (Elford, 1959). 


FUTURE WORK 

The application of the analysis proposed by Groves to the winds measured 
by the radio meteor method has been shown here to provide a powerful tool in 
the investigation of periodic atmospheric motions for the height range from 80 
to 100km. The demonstrated domination of this wind regime by solar influ- 
ence is only one of the many possible products of the analysis. It is hoped 
that attention can be given some time in the future to the determination of, for 
example, the magnitude of the lunar influence on winds at this height. For the 
present, attention is being given to the more detailed interpretation of the 
results already obtained. 
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figure 3a. Wind Energy Spectrum at 83km. February, 1961. 
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APPENDIX I 

Listing of the IBM 7094 FORTRAN IV 
Spectrum Analysis Program. 
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non 


NOR 


ENERGY SPECTRUM TO PLOTTER TAPE 


$ l LTGY b 

1IU Ju2 I 
SPAuot 
ILAtCulE IB JOB 

$IbJOU GU * LOG I C * MAP 

$ I B F T C ERG M94.XR7.LIST 

C ERG (ELFORD. ROPER* GROVES) METEOR WIND SPECTRUM ANALYSIS. 

C 

C ROUTINE SPECTRUM ANALYSIS. 

C OPTIMIZED DATA HANDLING FORTRAN IV PROGRAM. 

C DATA OUTPUT FROM 1 190A To BE ON TAPE 2. 

C MAXIMUM OF 2000 ECHOES* ALL FROM ONE MONTH. 

C 

c data on tape 2 consists of 

C HEADER card , REPRODUCED AI the top of each page of output. 

C FORMAT 12A6 

C OUTPUT TAPt DESIGNATION. FORMAT 6A6. 

C HEIGHT RANGE. ZMIN-ZMAX. FORMAT 1X2F4.0 

C NEXT CARD NORMALLY HAS A ZERO OR BLANK IN COLUMN 1. IF SEVERAL 

C SEGMENlS OF A SPECTRUM ARE REQUIRED. A 1 IS PUNCHED IN COLUMN 1. 

C STARTING FREQUENCY FOR PFRIODOGRAM. IN CYCLES/DAY (NORMALLY 0.5) » 

C END FREQUENCY (NORMALLY 4.0). AND THE FREQUENCY INCREMENT 

C (NORMALLY 0.05). FORMAT I1.3F5.2 

C DATA OUTPUT 119UA, MAXIMUM OF 2000 CARDS. 

C A BLANK CARD TO FLAG. END OF ECHO DATA. 

C IF A 1 HAa BEEN PUNCHED IN COLUMN 1 OF THE FREQUENCY RANGE 

C AND INCREMENT CARD. THEN FURTHER CARDS SPECIFYING THE SPECTRUM 

C INIERVALS AND INCREMENTS FOLLOW. EACH WITH A 1 IN COLUMN 1. 

c the last of these cards (if anyj must contain a zero or blank 

c IN COLUMN 1 TO TRIGGER THE OPERATOR REMOVE TAPE MESSAGE. 

C 

C ANY COMPLETE SET OF CARDS ABOVE MAY BE FOLLOWED BY ANOTHER 

C SET. COMMENCING FROM. AND INCLUDING. HEADER CARD. 

C 

C 

DIMENSION KHTI20) 

DIMENSION RESULT ( 12 ) 

DIMENsIONUl24.24) »A)24. 4e)»R(24.24)*P(24)*D(24).AC<24) 

DIMENSION s I GMA ( 24 ) 

DIMENSION LTHMM ( 2000 ) .EL3M( 2 000 ) .EM3M ( 2000 ) * ZM ( 2000) .VELMI 2000) * 

1 JOMI2000) 

DIMENSION T APE ( 6 ) 

DIMENSION FINIS) 126) 

COMMON A* R» NOP.ZMIN.MIN.ZMAX.MAX. SUM.NP .NQ.NR .NA0.NB0.NC0. 

1NA.NB.NC* AC. RESULT. PERIOD 

COMMON KHT.NHT.MY.MO 

This SECTION OF THE PROGRAMME READS PROCESSING PARAMETERS. 

60 uu CONTINUE 

iprint.u 

LOAD—1 

M»0 

ZERO.U.O 

NPASS.U 

NOGO.O 

NOP-O 

READ (2*31) RESULT 

31 FORMAT) 12A6) 

READ 12.31) TAPE 

NP« 1 
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NQ= 1 
NR = 1 

READ ( 2 *2 UMIN.ZMAX 

2 FORMAT ( 1X2FA.0 ) 

M I N=ZM I N 
MAX=ZMAX 

I F ( Z MAX -2 MIN-31 .0 ) 39 * 39 *40 
A v, .vRITE (3*41) 

a 1 FORMAT ( IX 3 7 HUE IGHT RANGE EXCEEDS ALLOWABLE 30 KM. /// 1X2 1HEXECUT ION 
1 I tRMINATEU. ///////) 

PRINT 41 
GOTO 30U 
39 KM 1(1) = MI N 

NhI = (/IMaX-AMIN)/2.0+1.0 

865 READ(2*866)IBIT*START*ENDIT*STEP 

866 FORMAT ( I 1 » 3F5 • 2 ) 

NUMPAS* ( END I T-START ) /STEP+l .6 
ENDIT = ENDIT-O.OU(J1 
CYCLE=START-6TFP 
LOAD=LOAD+ 1 
IF(L0AD)700*700>1 
7uu D070 I = 2 ♦ NH T 

KHT ( I ) =KHl ( 1-11+2 
7'-' CONTINUE 
NAU = 2 
NA=2 

NAHOLD=NA 
NB0 = 3 
NB= 3 
NCU = U 
NC = 0 
N = 24 
NE = 25 
N2 = 48 
NAUE = 3 
NAOT =4 
NAE = 3 
NA2 = 4 
NBUE=4 
NBOT = 6 
NBE = 4 
NB2 =6 
NCOE=l 
NCUT=0 
NCE= 1 
NC2 = 0 

PRINT 2000 * T APE 

20UU FORMAT l 1X4HL0AD 6A6*32H0N B6* AND PUSH START WHEN READY /1H1/1X) 
PAU^E 
REWIND 46 
CALL DNSHI 

1 WRITE(16)RESULT*KHT(1)*NHT 
PRINT 7193 *NUMPAS 

7193 FORMAT l 1X30HNUMBER OF PASSES THIS SEGMENT 16///////) 

I F ( LOAD ) 3 * 3 *999 

3 READ (2*4)UR.MP.MQ*J0.LTiMH*LTIMM»EL.EM,EL3*EM3*Z*LEVEL.LT. 

1NFL»NFM»VEL*MCS.MX 

4 FORMAT (F7.0.4 13* 12 *4F5.2»F4.0 *213 * 2 I 5 * F5 .0* 2 I 3 ) 

I F ( UR ) 999 * 999 ♦ 5 

5 I F ( Z— ZMAX ) 6 *6 • 3 
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n rs n o 


6 IK ( Z-ZMJN ) 3 ,7.7 

7 M = M + 1 
MY=MP 
MO = MO 
JOM ( M | = JO 

LTHMMIM ) =LT IMH*10U + LT I MM 

EL3M(M ) =EL3 

EM3MIM) =EM3 

ZM( M ) -L 

VELM(M) =VEL 

GOlO 3 

999 CYCLE=CYCLt+bTEP 
PER IOIJB24. O/CYCLE 
NA=NAHULD 
DO 3 8 I = 1 .2A 
P ( I ) =u.u 
AC I I )=0.U 
D038L=1 * 2 A 
Q(l'»L) = U.U 
38 CONTINUE 
SUM 1 =U , U 
SUM2=U.U 
SUM3=U.O 

NEXT CUMEb PROCESSING OF ECHO DATA TO PRODUCE COLUMNS D AND P» 
AND MATR IX U. 

I RUN=U 

6 u I RUN = I RUN+ 1 
J0= JOM ( IRUN ) 

L T HM = L I HMM ( IRUN) 

EL3=EL3M( IRUN) 

EM3=EM3M( IRUN ) 

Z = Z M ( IRUN ) 

VEL = VELM< IRUN ) 

CALL TIME IMY.MO.JO.LTHM.PERIOD.T ) 

DCL=EL3 

DCM=EM3 

DCN=SQR T l 1 .U-EL3**2-EM3#*2 ) 
b = ( 2 .0*Z-ZMAX-ZM IN ) / ( ZMAX-ZM I N ) 

S =S + U , UOUOG 1 
VEL=-VEL 
SUMP = U 
SUMO=0 
SUMR =U 
NCOUN T = 0 
SUMSAU = 1 

IFINAuT ) 11U.110.8A 
8 A D08K=2 .NAOT .2 

SUMSAO = SUMS AO + S**K. 

8 CONTINUE 
SUMS A =1 
D09K = 2 * NA2 * 2 

9 SUMSA=SUMSA+S**K 
SUMP = SUMP+SUMSA 

110 SUMP-SUMP+SUMSAO 
SUMSB0=1 

IFINBOT ) 1 3U . 1 30 . 1 1A 

1 1 A D011K*2»NBOT.2 

SUMSBO»SUMSBO+S#*K 

11 continue 
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n r\ n 


bUMi>B=l 
D012K = 2 »NB2 .2 
12 6UM5B=SUMSB+S*#K 
bUMQ=SUMQ+6UMSB 
13 o 6UMQ=SUMQ+SUMSB0 
SUMSC0=1 

IF1NC0T ) 1 60 » 160 * 1 35 
135 D014K=2 tNCOT.2 

6UMi>C0 = 6UMSCO + b*#K 

14 CONTINUE 
6UM6C=1 

D01 5K = 2 *NC2 »2 

15 6UMbC=oUMbC+b**K 
6UMR = SUMR + SUMSC 

16o SUMR = SUMR + bUMSCO 

WF=1.0/ ( ( DCL**2 )*SUMP+( DCM**2 )*SUMQ+< DCN*«2 )*SUMR ) 

SUM3=SUM3+WF*VEL**2 

D020K=1 *NAOE 

ncount=ncount+i 

20 D(NC0UNT)=DCL*(5**(K-1) ) 

D021K= 1 *NAE 

NCOUNT = NCOUNT + 1 

21 D(NC0UNT)=DCL*(5**(K.-1> ) * S I N IT) 

D022K= 1 »NAE 

NCOUNT = NCOONT + 1 

22 D(NC0UNT)=DCL*IS**(K-1) )«COS (T) 

D023K= 1 *NBOE 
NC0UNT=NC0UNT+1 

23 DINCOUNT )=DCM*(S#*(K-1) ) 

D024K = 1 »NBE 
NC0UNT=NC0UNT+1 

24 0 ( NCOUNT ) »OCM* ( S**( K-l ) ) N (T) 

D02 5K= 1 *NbE 
NC0UNT=NC0UNT+1 

25 D ( NCOUNT ) = DCM# I S*#( K-l ) ) #COS (T) 

D026<= 1 iNCOE 

NCOUNT = NC0UNT+1 

26 DINCOUNT > =DCN* ( S»* ( K-l ) ) 

D027K=1 *NCE 

NCOUNT sNCOUNT+1 

27 D ( NCOUNT ) =DCN»( S*#( K-l )) *SIN (T) 

D028K=1 *NCE 
NC0UNT=NC0UNT+1 

28 UINC0UNT)=DCN*(6*#(IC-1 ) I *C06 (T) 

D029 J= 1 *N 

P( JI=P< J)+WF*VEL*DU) 

29 CONTINUE 
D030J= 1 *N 
D030K* 1 *N 

Q( J.K)=Q( JtM+WF«D< J)*DU| 

30 CONTINUE 

I F ( IRUN-M160. 5894. 5894 

INVERSION OF Q. AND FORMATION OF COEFFICIENT COLUMN AC. 

5894 D0101 J* 1 *N 
D0101K=1 *N 
A( J*KI»Q(J*M 
101 CONTINUE 
D01 U2 J= 1 * N 
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n n n 


D01 02K=NE *N2 
I F ( J-K+N ) 108*107*108 
1»7 At J. K)=1.0 
GO TO 1 o2 
N8 A ( J *IC ) =0. 

102 CONt IMOE 

CALL MAlBlN (N.MISS) 

IF (MIG0I1103. 1103*1102 
11^2 CALL PAGE ( RESUL T . PE R I OD ) 

WRITE ( 3 * 1 104 ) N *N 

1104 FORMAT! 1H0 / 1X53H «#*** ERROR IN INPUT DATA HAS RESULTE 

ID IN MATRIX Qt I 3 * 1H * I 3 » 34H ) BEING UNSUITABLE FOR INVERSION.///// 
2 1 A 30HPROGRAMME CANNOT BE CONTINUED.///////) 

PRINT 1104 
GOT 0300 
1K3 CONTINUE 

D0103K= ] * N 
DO 1 03 J= 1 » N 

103 AC(K)=AC(<)+P( J)*Rl J*K) 

001040= 1 *N 

00 1 04K = 1 * N 

iUMl=BOMl +AC( J)*ACIM»0( J.K) 

104 CONTINUE 
D01 05 J= 1 *N 

SUM2=SUM2+AC( J)*P( J) 
lv5 CONTINUE 

SUM* ISUMl-l 2.0*SUM2 ) +SUM3 ) /FLOAT ( M-N ) 

D01 06 J= 1 * N 

1^6 SIGMAIJ) =SOR T (R(J»J)*SUM) 

PRELIMINARY OUTPUT. 


NOGO*NOGO+ 1 
IF I NOGO-1 171.71*72 
71 CALL PAGE (RESULT *ZtRO) 


WRITE (3*200) 

2uu FORMAT! 1H ////1X48H VARIATION OF UPPER ATMOSPHERE WINDS WITH HEIGH 
1I//1X51H BAoED ON GROVES ANALYSIS. WITH ERROR DETERMINATION////) 


WRITE 

201 FORMAT! 
1 BER OF 

2 1X27H 

3 1X15H 
WRITE ( 

2u2 FORMAT! 


13*201 ) M* N .MAX .MIN 

1H ////1X32H NUMBER OF METEORS PROCESSED 
INPUT PARAMETERS a 14//// 

HEIGHT RANGE* MAXIMUM I5»1X11H MIN 
KILOMETRES. ///> 

3*202) NAO » NA 

1X20H EAST - WEST PROFILE 2 I 3/ IX ) 


203 

204 
72 

206 


WRITE 

FORMAT ( 1X20H 


WRITE 
FORMAT! 1X20H 
CALL PAGE 


( 3*203 INBO.NB 
NORTH-SOUTH PROF I LE2 
I 3 * 204 ) NCO . NC 
VERTICAL PARAMETERS2; 
(RESULT. PERIOD) 


3/ 1 X ) 
3) 


WRITE (3*206) 

FORMAT! 1X2SH COLUMN MATRIX AC<X) /IX) 


■I 5////1X33H NUM 
MUM I 5 . 


U«0 


D0205 I * 1 . N 


JaJ+1 

I F ( J-501208.208.207 

207 JaO 

CALL PAGE (RESULT. PERIOD) 

WRITE (3*206) 

208 WRITE 0*209 ) AC ( I ) .SIGMAI I ) 

209 FORMAT ( 1X2SXF7.2*3XF7*1 ) 
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205 CONTINUE 

CALL PAGE (RESULT. PERIOD) 

CALL VARY 
NPASS=NPA6S+1 

WRITE (3. 299 ) NPASS 

299 FORMAT ( IX 12H END OF PASSI5) 

IPRINT= IPRINT+1 
I F ( I PR I NT-10 1281*280. 280 

280 IPRINTsO 
PRINT 299 .NPASS 

281 IF(CYCLE-ENDIT >999*300*300 
30U NhT6=6*NhT+6 

003011=1. NHT6 
3ol F I N I S ( I ) =0.0 

WRITE (16) (FINIS(l) *I*1»NHT6) 

I F ( I B I T >332.332.865 

332 END FILE 46 
REWIND 46 

PRINT 333* ( TAPEl IOUT | . IOUT= 1 .6) 

333 FORMAI ( 1X6HREM0VE/1X/1X 6A6.39HFROM B6 * AND HOLD AS FILE TA 

1PE. / // 1*16HTHEN PUSH START.//// IX. 10X12HMANY THANKS.////////) 

PAUSE 
GOTO 6000 
END 
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S' I B F T C VARY I- - 9 A ♦ XR7 * L I S T 
SUBROU I I N L‘ VARY 

C CALCULATES THE AMPLITUDE AND PHASE OF THE PERIODIC COMPONENT. 

C TOGETHER WITH THE MOST PROBABLE ERROR IN EACH. 

C PRINTS THEbE OUT AT 2KM INTERVALS OVER THE HEIGHT RANGE SPECIFIED. 

C 

DIMENSION AI2A.A8).R<2A.24) 

DIMENSION AC(2A) »AU(20) »PH(20) * ERAMP t 20 ) » ERPH ( 20 ) .RESULT ( 12 ) 
DIMENSION UO(20). EO(20). K.HTI20) 

COMMON A. R * NOP.ZMIN.MIN.ZMAX.MAX . SUM »NP . NQ . NR ♦ N AO » NBO .NCO . 

INA.Nb.NC. AC .RESULT .PER I0D 

Common mht .nht .my. mo 

N A o E = N A 0 + 1 
f.B0E=NB0+l 
NCUE =NCo+ 1 
KEND=u 
NS I GN = ~ 1 
I DENT = 1 

WRITE (3*597) 

597 FORMAT ( 1X39HEAST - WEST COMPONENTS OF THE MEAN WIND//) 

98 WRITE (3*87) (KHT( I )»I*1*NHT) 

87 FORMAT! 1X8HHEIGHT 15.1617) 

WRITE (3*871) 

871 FORMAT I IX ) 

I T = 0 
KA=KEND 

001 28KZ =M I N .MAX » 2 
I T= I T + l 
UO ( I T ) = 0, 0 

Z*Kl 

S=( 2.0*Z-ZMAX-ZMIN) / IZMAX-ZMIN) 

S=S+0. 000001 
D0106X= 1 ♦ NAOE 
KUA = K.+KA 

106 UO (IT) «UO (IT) +ACIKUA)*S**(K-1 ) 

S I GUO=0 .0 
D01 07K = 1 * NAOE 
D01 07L= 1 * NAOE 

1 u7 6 I GUO = R I GOO+S** ( K.-1 ) *S** ( L-l )«R(K. *L )*SUM 
EO ( IT) = SQRT ( S I GUO ) 

kend=ka+naoe 
numna=o 
si =0.0 

C0*0.0 
S I GCOS = 0 • 0 
SIGbIN=0.0 
SIGSC=0.0 

numna=numna+na 

NUSINsNP+NUMNA 

kstart»kend 

naend=na+i 

kend=kstart+naend 

doiiok*i *naend 

ks«k+kstart 

si»si+acus)#s#*(k-1) 

KC*KS+NUSIN 

no co«co+ac(kc)#s»*u-1) 

S I NSQJ=S I **2 
COSSQJ»CO**2 
SUMSO J> 5 1 N&OX-COSSQJ 
AU(IT)-SQRT (SUMSOJ) 
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F J = PER I OD 

IF (CO I 114*111 »1 15 

111 I F ( b 1 1113.113.112 

112 PHI IT ) =FJ/4.0 
GOTO 116 

113 PHI I T ) =F J*U. 75 
G0T0116 

114 PHI IT)=FJ*0.5+(FJ/6.28318)*ATAN (bl/COI 
GOTO 116 

115 PHI lT)=(FJ/fe.28318)*ATAN (bl/CO) 

116 IFIPHi IT ) >117.118.118 

117 PHI I T ) = F J + PH ( IT) 

118 DO 1 19K= 1 * NAEND 
DO 1 1 9L= 1 .NAEND 
Kb=K+KbT ART 

Lb = L+<bT AR I 
K.C = Kb + NUb I N 
LC=Lb+NUb I N 

bIGblN =bIGbIN +b*» I K-l ) *S**( L-l ) *R I Kb.Lb ) 
bIGSC =bIGbC +b**(K-l )*S**(L-1 )*R(KS.LC) 

119 bIGCOb =b I CiCOb +b**(K-l)*S**(L-l)*R(KC.L.C> 

PROD=2 • 0*b I »CO *b I GbC 


bIGPH = (CObbQJ*SIGbIN +b I NbQJ*S I GCOS 
bIGAMP =(bINbQJ*bIGbIN +CObbQJ*S I GCOb 
ERPHI IT ) = bQR T I bIGPH) *FJ/6. 28318 
121 ERAMPI I T ) =bQRT (bIGAMP) 

128 CONTINUE 

<end=kend+nusin 

WRITE 

FORMAT 1 1X6HME AN 
wRITE 

FORMAT l 1 X6HERROR 
WR I TE 
WRlTF 

FORMAT I 1X6HAMP 
WRITE 
WRITE 
WRITE 


-PROD )*bUM/bUMbQJ**2 
+PROD)-bUM/bUMbOJ 


123 


1 2 3u 


124 


125 


12 50 


872 


(3.123) (UOI IT),IT = 1* NHT ) 

16F7.0 ) 

(3.1230 ) I EO I IT). i T = 1 .NHT ) 

16F7.0 ) 

I 3.871 ) 

(3.124) (AU( IT).IT = 1.NHT) 

16F7.0 ) 

(3.1230) 1 ERAMPI IT ) » I T = 1 » NhT ) 

( 3.871 ) 

(3.125) ( PH ( IT ) ,IT = 1 .NHT) 

FORMAT ( lX5HPHAbF16F7. 1 ) 

WRITE ( 3 » 1250 ) ( ERPHI IT ) . I T = 1 .NHT ) 

FORMAT ( 1X5HERROR16F7. 1 ) 

WRITE (3.872) 

FORMAT l 1X/1X) 

WR I TE I 16) I DENT .MY.MO.PEHIOD.XHT ( 1 I .NHT . ( UO ( IT ) . I T = 1 .NHT ) 

1 l EOl I T ) . I 1 = 1 ,NhT ) . I A o I I T ) . I T = 1 >NhT ) . ( ERAMP ( I T ) . I T = 1 , NhT ) . 

2 ( PHI I T ) . 1 T = 1 .NHT ) • ( ERPHI I T ) . I T = 1 .NhT ) 


IFINbIGN) 129.131.133 
129 NAOE=NBOE 
N A = NB 
N b I GN = 0 
I DENT = 2 

WRITE (3.598) 

598 FORMAT ( 1 X 39HNOR TH-bOuT H COMPONENT b OF THE MEAN wIND/Z) 
GOT 098 

131 NAOE=NCOE 
NA = NC 
Nb I GN= 1 
I DENT = 3 

WRITE (3.599) 

599 FORMAT ( 1X39HVERT ICAL COMPONENT b OF THE MEAN WIND//) 
GO T 098 

133 RETURN 
END 
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n n r\ n r\ n 


SIHFTc MAlblN M94 » XR7 » L I i>T 

bORROuT [NF MA T b I N (N.Mlbb) 

i NV FRb I ON OF MATRIX OF ORDER N. UP TO 70 X 70. 

PROCEEOb VIA A METHOD OF GAUbb I AN ELIMINATION. DEbTROYING 
THE AUGMENTED MATRIX A IN THE PROCEbb. 

DIMENSION A (24.48)* XI24.24) 

COMMON A. X 
M I b b = — 1 
MM=2*N 
DO 15 I =2 *N 
70 11=1-1 

7 DO 15 J = 1 * I I 

8 IF ( A ( I .J> >9*15*9 

9 IF (ABO (A(J.JI)-ABb < A ( I . J ) I ) 1 1 * 1 0 . 1 0 

10 R = A( I » J ) / A ( J » J ) 

GO TO 130 

11 R = AU.J)/Al I .J) 

DO 12 K = 1 * MM 
B=AI J.< ) 

A ( J * K > = A < I *K) 

12 At I *M=B 
130 JJ=J+1 

13 DO 14 K. = J J * MM 

14 At I *K)=A( I *K.)-R*A(J*K) 

15 CONTINUE 

IF lABb ( A(N.N))-1.0E-10) 16*16*17 

16 M I bb= 1 
GOTO 29 

17 D028 J= 1 *N 
KK. = N + J 

X ( N.J ) =A ( N.K.K, ) /A ( N*N ) 

D028 I =2 » N 
J J = N- 1 +1 
L) = U. 

I I =N- 1 +2 
DO 2 5 iC = I I »N 
25 B=B+ A ( JJ*K.)*X(K»J) 

IF (ABb (At JJ.JJII-l.OE -10)16*16*28 
28 X ( JJ, J ) =( At JJ.KK ) -B) /At JJ. JJ ) 

29 RETURN 
END 
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noon 


SIBFTC PAGE M94.XR7.LIST 

subroutine page (result. period) 

T URN5 PAGE. NUMBERS IT. AND WRITES HEADING AS APPEARING 
ON RESULT CARD. 

DIMENSION RESULT( 12) 

DIMENSION A ( 24 .48 ) ♦ BI24.24) 

COMMON A. B. NOP 
NOP = NOP + 1 
I F ( PER I OD ) 2 »2 . 3 

2 WRITE ( 3» b (RESULT .NOP 

5 FORMAT ( 1H1/ 1X12A6. 1uA17HoPEC TRUM ANAL Yi> I 53A4HPAGE I 5/ / / ) 
GOT04 

3 FREQ=2A.O/PERIOD 

WRITE (3.1 (RESULT .PERIOD .NOP. FREQ 

1 FORMAT ( 1H1/ IX] 2A6* 1UX6HPER JODF8.2 > 5HHOURS 7X4HP AGE I 5// 

1 1X82X9HFREQuENCYF8.A. 12H CYCLES/DAY./// ) 

4 RETURN 
END 
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S I BFT C TIME M9A»XR7*LIi>T 

subroutine time (MYtMo. jo. lthm, period. t i 
calculates time of echo with respect to input periodicity. 

LT I MH = L T HM / 1U0 
LT I MM = L T HM-L T IMH*1U0 

TMINI T = ( JO-1 )*1440+LT1MH*60+LT IMM-8 
ThOUR=TMINI T/60.0 
NEWDAY=ThOUR/PER iod 

T=( THOUR/ PERIOD-FLOAT ( NT *DA Y ) 1*6.28318 

return 
END 
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$ I 8MAP 

DNSHl 

H94/2 *XR7 


ENTRY 

DNSH I 

DNSH I 

SDHB 

6 


tra 

L *4 


END 
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APPENDIX IA 

Listing of the IBM 7094 FORTRAN IV 
Master Tape Production Program. 
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$ 1 BFT C COPY M94.XR7.LIST 
SUBROUTINE COPY 

D I MENS I ONUO ( 20 ) . EO ( 20 ) . AU < 20 ) . ERAMP ( 20 > .PH! 20 > . ERPH! 20 ) .RESULT! 12 ) 
REWIND 47 

READ! 17 IREbULT .KHT .NHT 
WRITE! 16 )REoULT.KhT.NhT 

1 READ! 1 7 ) 1UENT . M Y » MO . PER I OD . K . L . <uO( I R ) . I R= 1 . NH T ) . 

1 l E 0 ( I R ) » I R = 1 » NHT ) 1 1 Aul I R ) » I R = 1 » NHT ) .( ERAMP ( I R ) » I R = 1 » NhT ) . 

2 ( PH! IR ) »IR=1 .NHT ) . ( ERPH ( I R ) . I R=1.NHT ) 
wR I TE I 16 ) I DENT .MY ,M 0 .PERIOD. K »L. ! UOl IR ). IR = 1 .NHT ) . 

1 ( EO! IR I . IR=1 .NHT) , ( AU( IR) . IR=1.NHT ) . ( ERAMP! IR ) . IR=1 ,NHT ) . 

2 (PH! IR ) . IR=1 »NHT ) , ( ERPH ( I R 1 . I R = 1 .NHT ) 

IF (MO) 2. 2*1 

2 REWIND 47 
RETURN 
END 
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S1BSYG 

SID J02 T RGR BLDG 1. COPY ERG OUTPUT TAPES ONTO MASTER TAPE. 

SPAUSE 

SEXECUTE I B JOB 

$ I B JOB GO. LOGIC. MAP 

SIBFTC MASTER M94.XR7.LIST 
C 

c copies output tapes from erg onto one file tapf. 
c 

DIMENSION tape! 14) 

DATA l I APE I I ) » I = 1 ♦ 14 ) /6HAMU850.6HAM2354.6HAM2283.6HAM2321 .6HAM2 2 24 
1 * 6H AM2 1 7 7 * 6HAM2 24 7 * 6HAM2 1 07 * 6HAM1 0 33 ♦ 6HAM2 266 * 6PIAM24 36 » 6HAM2 1 5 3 ♦ 

1 6HAM22 9 9 . 6 HAM 2 391/ 

CALL DNSH1 
PRINT 1 

1 FORMAT ( 1X31 HLOAD I APE AM286u ON Ufc» AND /IX) 

D03 I = 1 » 14 

PRINT 2 . TAPE! I ) 

2 FORMAT ( IX 1UHL0AD TAPE A6.23H ON B7* THEN PUSH ST AR 1 /////// ) 

PAUSE 

CALL COPY 

3 CONTINUE 
KH T sU 

Nh I =0 

w R I f F < 16) I I APE I I ) . I = 1 »12 ) .KhT .NhT 
LND Fill. 46 
REMIND 46 
PRINT 4 

4 FORMA I ( 1X61HREM0VE AM2850 FROM b6 » AM2391 FROM b7* AND Push o I AR I 
1 TO FX I I /////// ) 

return 



S1BMAP DNoHl M94/? * XR7 *L 1 ST 

DNSHI ^AVf. ( 1 *2»3 »<m5»6W) I 

SDMB 6 

^DHB 7 

l HA 1*4 

END 
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APPENDIX H 

Listing of the IBM 7094 FORTRAN IV 
Plotter Program. 


46 


3- I RSYS 

SHI JU2T HGR BLDG 1. WRUEb PLOTTER T APES FROM ERG MASTER. 

$PAUSE 

3 E A L C U I L IB JOB 

1 IB JOB GO* LOGIC. MAP 

$ I BF T C PLul M94.XR7.L1GT 

C RL J I R"olI,.E i -IK. IV METEOR WIND ENFRGY PER IODOGRAM. 

C I NCLuJLo I H I ERPOLA T I ON ROUTINE. 

C READS ERG MASTER TAPE (B7) 

C >vi< I I FS PLOTTER TAPE ( B6 ) 

C II,.;. PRuVIoION FOR REPORT SUE OUTPUT GRAPHS. 

C READS TAPE 2 DATA AS FOLLOWS 

C YEAR. MONTH FOR WHICH OUTPoT |o DESIRED. FORMAT 213 

C MODE » PUj I I I ON oPEC IF ICA T ION.HE IGHT FOR WHICH SPECTRUM IS REQUIRED. 

C jRUINAlF NORMAL 1 / A 1 ION FORMAT 2I1.I3.I6 

C MODE S L L E C i 1 0 N 

C Total E hER GY . Ew*»2 + Ns»*2 + VERUCAL**2 

C 1 Ew AMPLIIUDl 

C 2 Ew PHAot 

C 3 Ns AMPLITUDE 

C 4 NS PHASE 

C 6 VERTICAL AMPLITUDE 

C 6 VERTICAL PHASE 

C 7 EW AMPL I T Jt)F»*2 

C b AMPL I 1 UDE **2 

C 9 vfRTiCAL AMPL I ToDE**2 

C oP lu o I X GRAPHS CAN BE PLOTTED AT REDUCED SCALE (OUTPUT 

C APPROPRIATELY POu I I 1 ONED ON PAGE 8 X 11). IF NO POSITION 

C uPECIKICATlON. OUTPol GRAPH WILL FILL FULL PAGE. 

C IF NO ORDluATE NORMALIZATION IS SPECIFIED. OUTPUT GRAPH WILL 

C AulOMA I ICALLY FILL FRAME. 

C CON T I Nut wITh AS MANY OF THESE CARDS AS REQUIRED 

C A BLANK CARD RETURNS CONTROL TO THE START OF THE PROGRAM 

C FOR MULTIPLE MONTH HANDLING. 

C A NEGATIVE HEIGHT FLAGS END OF DATA. AND OPERATOR REMOVE TAPE 

C MESSAGES. 

DIMENSION T APF ( 12 ) >X ( 2000 ) .Y12000) 

COMMON/ AMPS/ UO ( 20 ) ♦ EO ( 2 U ) . AU ( 20 ) »ERAMP (20).PH(?0).FRPH(20>. 

1RESULT ( 12 ) 

COMMON/ XY/NJ .NX ( 200U ) .NY ( 2000 ) 

C 1 = 0. 1 
C=C 1 

T WOP I =6 .2832 
PRINT 1 

1 FORMAT ( 1X21HLOAD A M 2316 ON B 6 / 1 X / 1 X2 2H AND A M 2850 ON p 7 / 1 

IX/ IXI 8 H 0 CRATCH TAPE ON B4 /I 

I X / 1 X2S H AND PUSH START WHEN READY///////) 

PAUSF 
REWIND 44 
REWIND 46 
REWIND 47 
CALL DNSHI 

2 I OUT = 0 

wR I TE ( 3 . 1000 ) 

10UU FORMAT ( 1 H 1 / 1 x ) 

C 

C DETERMINE MONTH TO BE PROCESSED 

C 

READ! 2 »4 IMYEAR .MONTH 
4 FORMAT (213) 

C 
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C SEARCH I APE ON B7 FOR DATA REQUIRED 

C 

CALL oLARLHlMYEAR *MON T H ) 

C 

C COPY REQUIRED DATA ONTO BA (SCRATCH TAPE) 

C 

CALL COPY 
C 

C READ PROCESSING PARAMETERS 

C 

3 READ ( 2 « A I MODE .XPUT * IH.NYMAX 
j u S L' = b 
NlYRC=2 
I K ( IH ) 300.2*5 
5 REMIND A A 
C 

C SET UP FRAME FOR OUTPUT GRAPH 

C 

CALL FRAME 
IFUPUT >71*71*70 
7- CALL SIX! kPUT ) 

71 CALL IbDRAwl 
C 

C READ AND PROCESS ERG OUTPUT 

C 

READ ( 1A ) RESULT * KHT *NHT 
I ( =0 

1 = 1 IH-XmT )/ 2+l 
7 I T= I T+ 1 
X ( I T ) sO.u 
Y ( I T ) =0.0 
IFIMODE ) 10*10. 11 
lw D09o J= 1*3 

READ( 1 A ) IDEN T , MY. MO .PERIOD »K*L. ( UO ( IR ) * I R= 1 *NHT ) . 

1 l E0( IR ) * IR=1 *NHT ) . ( Au( IR ) . IR=1 *NHT ) . ( ERAMP (IR)*IR=1 *NhT). 
2 ( PHI IR ) . I R = 1 *NHT ) . ( ERPH ( I R ) * I R = 1 * NHT ) 

IFl MO )1U0. 100*80 
8u Y ( I T ) = Y I I T ) +AU ( I ) **2 
9 u CONTINUE 

X( I T )=2A.Q/PERI0D 
G0T07 

11 GOTO! 12*13*1A*15*16* 1 7 * 1 8 . 1 9 * 20 ) .MODE 

12 JUSE = 1 
NT YPE=0 
GOT 02 1 

13 JUSE= 1 
NTYPE=1 
G0T021 

1A JUSE=2 

ntype=o 

G0T021 

15 JUSE=2 
NTYPE=1 
G0T021 

16 JUSE=3 
NT YPE = 0 
GOT021 

17 JUSE=3 
NTYPE=1 
G0T021 
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r> n 


18 Joot=l 
NTYPL=2 
Gu ( 02 1 

19 JujL=2 
N I YPE=2 
GUT02 1 

2u JUSE=3 
N I YPF=2 

21 D09 J = 1 < ' 

READ! 1 A ) IDLNT.MY.MO*PERIOD*K*L*(uO( IK ) * IR=1 *NHT ) * 

1 l E01 I R I * IK = 1 »NhT ) * ( AUl IR ) * I R = 1 *NHT ) . ( ERAMP ( ! R ) » I R = 1 * NhT ) » 
? ( PH< IK ) . IR=1 *NhT ) * ( ERPH ( IR ) ♦ IR = 1 .NHT ) 

I E ( MU ) 1 0 u » 1 u o • 8 

8 IK J-JUOE >9*22*9 

22 IFINTYPE-l >23*24*25 

23 Y 1 I T ) =AU( I ) 

G0TU9 

24 PHASE= ( PH< I) /PERIOD) *TWOPI 
Y ( I T ) =PHASE 

G 0 T 0 9 

25 Y ( I T ) = AU ( I >**2 

9 CON 1 I NU F 

X ( I r ) = 24 • U/PER IOD 
GO T07 
C 

C bEf UP AKRmY FOR PROCESS I NO BY IoDRAtv. 

C LINtAR INTERPOLATION BETWEEN DATA POINTS 

Ob)ED TO PROVIDE FILLERS TO OVERCOME PLOTTER PEN INFRTIA. 

additional filler routinfs included in ipdraw. 
i^ n = i T-i 

B IGY=0. 0 
DO 1 02 J = 1 * IT 

I F ( Y ( J ) -R IGY ) 1 U2 * 1U2 * 1 0 1 

101 B I G Y = Y ( J ) 

1 02 CONTINUE 
I B I GY=b I G Y+ 1 . 05 
I F ( I B I GY -NY MAX 1202*202*201 

2 0 1 NYMAX= I B I GY 

202 YMAX=NYMAA 
XMAXrX l IT ) 

XMULT=40oU. 0/XMAX 
YMUL T = 4000. 0/Y MAX 
NX ( 1 ) =X ( 1 ) *XMULT 
NY ( 1 ) =Y ( 1 ) *YMt)LT 
N I = 1 
NJ = 1 

1 o 5 N I =N I + 1 

IFINI-IT ) 1U9»109*107 
lu9 NJ = NJ + 1 

NS I GNX= ( ABS ( X ( N I ) - X ( N I — 1 ) ) ) / ( X ( N I ) — X(NI-1))+C 
YGR AD= l Y ( N I (-Y(NI-l) ) / I X ( NI ) — X ( N I — 1 ) )*( YMULT/XMULT ) 

NPO I NT = ( XKULT* < X(NI)-X(NI-l)))/2.0-1.0 + Cl 
D01 06 1 = 1* NPO I NT 
NX ( NJ ) =NX ( NJ-1 ) +2*NSIGNX 
Ybl T=NY (NJ-1 ) 

NY(NJ)=YBIT+YGRAD«2.0+C1 
IFINJ-2000) 108*107*107 
108 N J=N J+ 1 
106 CONTINUE 

NX ( N J ) = X ( N I )*XMULT 
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I 


NY ( NJ ( = Y ( N I )*YMULT 
I F ( NJ-2 OuO ) 105*107.107 
K7 CONTINUE 

MAXX=XMAX+C 
MAXY=YMAX+C 
IFIKPUT 173*73.72 
7? CALL S I X ( aPUT ) 

73 CALL ItJDRAW 
I OoT = I OUT + 1 

I F ( IOuT -4)1071 * 1071 * 1070 

107 0 lOUTrO 

wR I TE ( 3 * 1 000) 

1071 CONTINUE 

wR I TE I 3* 104 )RE OOLT . I H *MAXX *MAX Y * I B I GY *KPUT * NI *NJ. JUi>E *N TYPE .NODE 
lo4 FORMAT ( 1X12A6*/1X/1X6HHEIGHT I 6 * / 1 X/ 1 X6HXMAX I 6 * / IX / 1 X6H YMAX '16/ 
11X/1X6HBIGY I6/1X/1X6HXPUT 16/ 

2 1X/1X2HNI 15. 3H NJ 15* 5H JUSE 15* 6H NTYPE 15. 5H MODE 15 //IX) 

GOT 03 

3uo REWIND 44 
REWIND 46 
REWIND 47 
PRINT 301 

3ol FORMAT ( 1A27HREMOVE A M 2316 FROM B6/ 1X/1X24HAN0 A M 2850 F 
1 ROM B7/ IX/ 1X22HAND PuJH i>TART TO EXIT///////) 

PAUJE 7 

RETURN 

END 
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SIBMAP 

I BOR A »» 

M A / 2 * X R 7 


IBDkAft 

oAVf 

t 1*2»3»A»5*6*7)I 


* 


A AIT COMPLLT ION 

OF ALL IOEX ACTIVITY 

CHAN 

LOU 

2 

NO OF CHANNELS TO CLEAR 

Two 

FQu 

2 



A A f 

I AO 9 A 

DO TWICE 

CLF AM 

AX I 

CHAN • 2 


NLa i 

AtC 

LOAD % 2 . 

GET CHANNEL INDEX TO IR1 


l l r* 

• C HX A C 

TEST CHANNEL 


I K A 

*-l 

BUSY 


T IX 

NEXT *2,1 

GO TO NEXT CHANNEL • 


I IX 

CLEAR »A»1 

ONCE MORE TO MAKE SURE 


flNB 

ZERO 

TURN TRAPS OFF 

* 


ALL I/O NOW UNDER 

CONTROL OF IbDRAW 


^DLB 

6 



AX I 

120UU » A 



bTZ 

REC+1200C *4 



1 I X 

*-1*4. 1 



LX A 

1 tl 



rxi 

*+l * 1 *-l 



AX T 

0.2 



AX T 

1 2 OuO . 4 



o 1 L 

COUNT 



CAL 

belect 



GL w 

REC 


GO 1 

*7BB 

6 



RCHb 

CMD 



I COB 

* 



I RCB 

* + 2 



1 R A 

I UBF 



Bi>RB 

6 



tcob 

* 



ira 

GUI 


tube 

nop 




CLA 

NX ,2 



ALb 

IB 



ORA 

NY $ 2 



oL W 

T IME 



ORA 

READ 



oL w 

REC 


GO 1 A 

WTBB 

6 



RCHB 

cmd 



tcob 

# 



I RCB 

* + 2 



tra 

OUTR 



BbRB 

6 



tcob 

* 



tra 

GUI A 


OUTR 

NOP 




AX T 

12U00.4 



bTZ 

REC+ 1 2000 *4 



T I X 

1 »4 ♦ 1 



AX T 

1 2 000 » 4 



CLA 

time 



ORA 

REDROP 



5LW 

REC 



T X I 

*4-1 * 4 * - 1 b 


GO 

6TZ 

bwi 



CLA 

N Y » 2 



bTO 

TIME2 



CLA 

N Y + 1 » 2 
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i>UB I T 


ADD IT 


OUT 


STO 

TIME 

bUB 

NY. 2 

bTO 

bSP 

FI 

bUB 

= 11 

TM I 

OUT A 

CLA 

FI 

TPL 

ADDI T 

CLA 

FI 

ADD 

= 11 

r pl 

OUT 

bTO 

FI 

TZE 

OUT 

CLA 

T I ME 2 

bUh 

= 11 

b T 0 

T I ME2 

CLA 

NX ,2 

ALb 

18 

ORA 

T I ME 2 

ORA 

REDROP 

bLW 

FOG 

bLW 

REC+12000.A 

T IX 

5UB I T . A » 5 

bX A 

* + 2.1 

TbX 

OUTK . 1 

AX T 

**. 1 

tra 

GUBI T 

CLA 

FI 

bUB 

= 11 

bTO 

FI 

TM I 

OUT 

UE 

OUT 

CLA 

T IME2 

ADD 

= 11 

bTO 

r IME2 

CLA 

NX ,2 

ALb 

18 

ORA 

T I ME 2 

ORA 

REDROP 

bLW 

FOG 

bL w 

REC+ 12000 ♦ A 

I IX 

ADDI T .A. 5 

oXA 

* + 2 * 1 

[ bX 

OUTK. 1 

AX T 

*# » 1 

tra 

ADDI T 

bTZ 

UW69 

CLA 

NX ,2 

ALb 

18 

ORA 

time 

ORA 

REDROP 

bLW 

FOG 

bLW 

REC+12000.A 

T I X 

* + A » A . 5 

bX A 

* + 2 . 1 

f bX 

OUTK. 1 

AX I 

** . 1 

CLA 

GW69 

tze 

* + 2 

tra 

NAT2 


STL Sw/69 

TRA OUT + 1 

OuTa CLA NX, 2 

ALb 18 

ORA TIME 

ORA REDROP 

bL w FOG 

bLw REC+12000,4 

T I X *+4*4*5 

bXA *+2*1 

IbX OUlK.l 

AX T **,1 

NAT2 T X I *+l,2*-l 

T I X GO . 1 , 1 

PX A 0,4 

bTO FI 

CLA =12u00 

bUB FI 

ALS 18 

oTD CMD1 

G02 WTBB 6 

RCHB CM01 

TCOB * 

TRCB *+2 

IRA END 

BbRB 6 

rcoB * 

TRA GO 2 

END WEFB 6 

* RE I Nb T AT E CONTROL OF I/O BY IOEX 

ENABLE NULL 

TCOA * 

RDCA 

TCOB * 

RDCB 

NZT . IRPSW 

ENB* .TRAPX 



RETURN 

I BOR AW 

OUTx 

NOP 


OUTK2 

WTBB 

6 


RCHB 

CMD 2 


I COB 

* 


IRCB 

* + 2 


IRA 

END1 


BbRB 

6 


TCOB 

* 


TRA 

OU TK2 

END1 

AX T 

120U0 *4 


bTZ 

REC + 12000 *4 


T IX 

*-l *4* 1 


AX T 

12000*4 


CLA 

FOG 


ORA 

REDROP 


bLW 

REC 


T X I 

* + 1 * 4 * - 1 5 


IRA 

ltl 


BCI 

1*0013)) 

FMT 

BCI 

1* (IX. ( 1 

CMD 1 

1 OR T 

REC + 1 0 * **♦ 

CMD 

I ORT 

REC* *5 

CMD2 

I OR T 

REC* * 12000 
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CMDlu 

XY 

1 

NX 

NY 

select 

REDROP 

UP 

DOWN 

READUP 

READ 

REC 

YAWL 

El 

SW 1 

TIME2 

SW69 

FOG 

COUNT 

TIME 


LOAD 

ZERO 


IORT REC * » 1 50 

CONTRL XY 

UbE X Y 

Bbb 1 

BSb 2000 

BSb 2000 

USE PREVIOUS 

OCT 600000040000 

OCT 600000600000 

OCT 000000100000 

OCT 000000200000 

OCT 600000500000 

OCT 600000400000 

BSS 12000 

BSb 100 

Bbb 100 

BSb 1 

BSS 1 

BSb 1 


Bbb 

10 




Bbb 

1 




Bbb 

1 




AXC 

4*1 

CHANNEL 

E 

INDEX 

AXC 

3*1 

CHANNEL 

D 

INDEX 

AXC 

2*1 

CHANNEL 

C 

INDEX 

AXC 

1*1 

CHANNEL 

B 

INDEX 

AXC 

U*1 

CHANNEL 

A 

INDEX 


NULL 

PZE 

END 


54 


non 


SIBFTC FRAME M94 *XR7 * L I i> T 

bUBROUT INE FRAME 

GETS UP FRAME FOR OUTPUT GRAPH 

COMMON/XY/NJ.NX (2000 ) .NY ( 2000 ) 

NJ=1632 

NX ( 1 ) =U 

NY ( 1 ) =4030 

D0500 1=2* 403 

NX ( I ) =NX ( 1-1 ) + 10 

NY ( I )=NY( 1-1 ) 

5Uu CONTINUE 

D05000I =403*408 
NX ( I + l )=NX< I ) 

NY ( I+l ) =NYl I ) 

50uU CONTINUE 

D0501 1 =409*811 
NX ( I ) =NX ( I — 1 ) 

NY ( I) =NY ( I - 1 ) - 1 0 
501 CONTINUE 

D05UU1 I =811*016 
NXl 1+1 ) =NX( I ) 

NY( I+l )=NY< I ) 

50ol CONTINUE 

005021=817*1218 
NXl I > = NX ( 1-1 ) - 1 0 
NY ( I ) =NY ( 1-1 ) 

5u2 CONTINUE 

D05002 1=1218.1223 
NX ( I + l ) =NX ( I ) 

NY ( I +1 ) =N Y ( I ) 

50U2 CONTINUE 

005031=1224.1626 
NX ( I ) =NX ( 1-1 ) 

NY ( I ) = N Y ( 1-1 1 + 10 
5 u 3 CONT INUE 

D050U3I=1626*1631 
NXl I + l )=NX( I ) 

NY ( I + l ) =NY ( I ) 

5uu3 CONTINUE 
RETURN 
ENO 
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n n o 


SIBFTC SEARCH M94.XR7.LIST 

SUBROUTINE SEARCHIMYEAR. MONTH) 

SEARCHES TAPE ON B7 FOR APPROPRIATE YEAR. MONTH. 

COMMON /AMPS/UOI 20 ) . EO ( 20 ) . AU ( 20 ) . ERAMPI 20 ) .PH ( 20 ) .ERPH ( 20 ) . 
1RESULT (12 ) 

luu READ(17)RESULT.KHT.NHT 

if(kht n.i.2 

1 REWIND 47 
G0T01U0 

2 READ ( 17 ) I DENT .MY .MO.PER IOD.K *L* ( UO ( I R I . I R* 1 .NHT ) . 

1 ( E0( IR > . I R=1 .NHT I . ( AU( I R ) . IR = 1 .NHT ) . ( ERAMP ( I R ) . I R= 1 . NHT ) . 

2 ( PH( IR ) . I R = 1 .NHT ) . ( ERPH ( IR ) . I R= 1 .NHT ) 

IF(MO) lOu. 100.3 

3 I F ( MY -MY FAR >2.4.2 

4 I F ( MO-MON T H ) 2 . 5 » 2 

5 BACKSPACE 47 
BACKSPACE 47 
RETURN 
END 
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non 


SIBFTC COPY M94.XR7.LIST 
iuBROOIlNE COPY 

COPIES DATA REQUIRED BY PLOT FROM B7 ONTO BA (SCRATCH TAPE) 

COMMON/AMPS/UOI 20 ) .EOI 20) ♦ AU ( 20 ) »ERAMP( 20 ) .PH I 20) »ERPH<20> * 
1RESULT ( 12 ) 

REWIND 44 

READ! 17 ) RESULT *KHT .NHT 
wRUE(14)RESULT.KHT.NHT 

1 READ! 17)1 DENT »MY »M0 »PER I OD .K * L » < UO ( I R ) * I R= 1 * NHT ) . 

1 ( EOI IR ) . I K= 1 .NHT ) . ( AU< IR) . I R= 1 »NHT ) , (ERAMPI IR),IR=1.NhT). 

2 ( PH I I R ) * I R = 1 *NHT ) . lERPUl IR) . IR=1 .NhT ) 

WRITE! 14) l DENT »MY .MO.PERIOD.K.L* (UOI IR) * I R* 1 *NHT ) . 

1 ( EOI IR I . IR=1 .NHT ) . ( AU ( I R ) . IR = 1 .NHT ) . ( ERAMP ( 1 R ) » I R = 1 * NHT ) . 

2 ( PHI IR ) . I R= 1 .NHT ) . (ERPHl I R ) » I R=1 .NHT ) 

I F I MO 12.2.1 

2 END FILE 44 
REWIND 44 
RETURN 
END 
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n nno 


JIBFTC SIX M94.XR7.Lt ST 
SUBROUTINE SIXUPUT) 

REDUCED SCALE GRAPHS IREPORT SIZE. PAGE 8 X ID 
SIX PER PLOTTER TABLE SHEET. 

COMMON/XY/NJ .NX ( 2000 ) .NY (2000) 

AX1=1280.0 
AX2*2560.0 
AY=1720.0 
XMULTs0.24 
YMULTsO* 32 

GOTO! 1 *2*3*4*5»6) .KPUT 

1 ADDX=U.O 
ADD Y = (J • 0 
GOT07 

2 ADDX=u.O 
ADDYsAY 
GO I 07 

3 ADDX=AX1 
ADDY=o.O 
GOT07 

4 ADDXsAX 1 
ADDY=AY 
GOT07 

5 ADDX=AX2 
ADDY =0.0 
GOT07 

6 ADDX=AX2 
ADDY=AY 

7 D08 I = 1 » N J 
R = NX ( I ) 

NX ( I )=R#XMULT+ADDX 
R = NY ( I ) 

NY ( I )=R*YMULT+ADDY 

8 CONTINUE 
return 

END 
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i I BMAP 
DNSHI 


DNbUI M94/2*XR7*LIST 

i>AVE ( 1 *2 * 3 *4 >5 *6 • 7 ) I 

DDHB 4 

^DHH 7 

fRA 1*4 

END 
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